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The fundamental theory of crystal twinning has been long established, leading to 
a significant advance in understanding the nature of this physical phenomenon. 
However, there remains a substantial gap between the elaborate theory and the 
practical determination of twinning elements. This paper proposes a direct and 
simple method - valid for any crystal structure and based on the minimum shear 
criterion - to calculate various twinning elements from the experimentally 
determined twinning plane for Type I twins or the twinning direction for Type II 
twins. Without additional efforts, it is generally applicable to identify and predict 
possible twinning modes occurring in a variety of crystalline solids. Therefore, 
the present method is a promising tool to characterize twinning elements, 
especially for those materials with complex crystal structure. 



1 . Introduction 

Crystal twins are commonly observed during solidification, 
deformation, solid-state phase transformation and recrys- 
tallization in a variety of crystalline solids with low stacking 
fault energy. Often, these features occur on the nanometre to 
micrometre scale, and they represent a particularly symmetric 
kind of grain boundary, giving rise to a much lower level of 
interfacial energy than general grain boundaries. As an 
underlying mechanism for microstructural changes, crystal 
twinning has acquired great importance in fields such as 
metallography, mineralogy, crystallography and physics. 

Early efforts to define crystal twins were based on the study 
of deformation twinning. By convention, a deformation twin is 
a region of a crystal that has undergone a homogeneous shape 
deformation (simple shear) in such a way that the resulting 
structure is identical to that of the parent (matrix), but 
differently oriented. A twinning mode is fully characterized by 
six elements: (1) K t - the twinning or composition plane that is 
the invariant (unrotated and undistorted) plane of the simple 
shear; (2) rj 1 - the twinning direction or the direction of shear 
lying in Ki, (3) K 2 - the reciprocal or conjugate twinning 
plane, the second undistorted but rotated plane of the simple 
shear; (4) r\ 2 - the reciprocal or conjugate twinning direction 
lying in K 2 \ (5) P - the plane of shear that is perpendicular to 
K 1 and K 2 and intersects K 1 and K 2 in the directions r] 1 and rj 2 , 
respectively; (6) y - the magnitude of shear. Moreover, the 
orientation relationship between two twin-related crystals can 
be specified by simple crystallographic operations: a reflection 
across K l or a 180° rotation about the direction normal to K x ; 
or a 180° rotation about r\ x or a reflection across the plane 
normal to r\\. According to the rationality of the Miller indices 
of Ki, K 2 , r\i and r\ 2 with respect to the parent lattice, crystal 



twins are usually classified into three categories: Type I twin 
(K 1 and rj 2 are rational), Type II twin (K 2 and r\ l are rational) 
and compound twin (K lf K 2 , r\ x and r\ 2 are all rational). 

The classical definition and description of deformation 
twinning have been further extended to describe other twin- 
ning processes associated with phase transformation and 
recrystallization. Notably, the concept of transformation 
twinning is widely adopted for the elucidation of structural 
changes during martensitic transformation. Although the 
formation of twinned martensitic variants is driven by a 
deformation from the parent phase and may not have any 
relation to the simple shear deformation defined by the 
twinning shear, the detwinning process can be well predicted 
by these elements, especially for the newly developed ferro- 
magnetic shape memory alloys (Gaitzsch et al, 2009; Wang et 
al. , 2006; Li et al. , 2010). In such a case, the twinned martensitic 
variants always form regular arrays of alternate lamellae with 
fixed thickness and the twin boundaries are highly glissile, 
where the detwinning shear determines the shape memory 
performance. 

For many years, constant attempts have been made to 
determine twinning elements of crystalline materials from the 
knowledge of crystal structure, because of their importance 
for insight into possible twinning modes and resultant orien- 
tation relationships of twinned crystals in the context of 
microstructural manipulation. A systematic theory was 
developed by Kiho (1954, 1958) and Jaswon & Dove (1956, 
1957, 1960) based on the minimum shear criterion, and later 
completed by Bilby & Crocker (1965) and Bevis & Crocker 
(1968, 1969). It provides the general expressions - valid for all 
crystal structures - to predict the twinning elements for both 
Type I and Type II twins with a known twinning shear. 
However, in a practical determination of unknown twins, it is 
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only feasible to resolve the possible twinning plane for 
Type I twins or the twinning direction rji for Type II twins by 
means of transmission electron microscopy (TEM) or scan- 
ning electron microscopy/electron backscatter diffraction 
(SEM/EBSD). In other words, with the given general 
expressions, one always suffers from insufficient information 
to derive the unknown twinning elements, especially the 
twinning shear. As a common practice, laborious geometrical 
examination of the lattice correspondence of the stacking 
planes parallel to the twinning plane has to be conducted. 
Such a process becomes particularly difficult when the twin- 
ning plane and the shear plane are irrational and the crystal 
structure is complicated. Hence, there exists a substantial gap 
between the elaborate theory and the practical determination. 

In this paper, we present a complete method to find all 
twinning elements for the three classical types of twins, based 
on the assumption that a simple minimum shear operation 
transforms the lattice points of a crystal into their counterpart 
twin positions. The initial inputs are simply the crystal struc- 
ture and the experimentally determined K 1 (Type I) or r] 1 
(Type II). As a general method applicable to any crystal 
structure, it may facilitate future characterization studies of 
crystal twinning. 

2. Methodology 

2.1. Determination of twinning mode 

For a twinned crystal, the crystallographic orientations of 
the twin and its parent can be experimentally determined with 
SEM/EBSD or TEM. In the case of SEM/EBSD examination, 
the orientation of a crystal with respect to the macroscopic 
sample coordinate system is usually characterized in terms of 
three Euler angles. The misorientation between the twin and 
the parent is then calculated from their Euler angles, and 
expressed by a set of rotation angles and the corresponding 
rotation axes (Cong et al., 2006, 2007). According to the 
definition of twin relationships mentioned above, there exists 
at least one 180° rotation. If the Miller indices of the plane 
normal to the 180° rotation axis are rational, the twinning 
mode belongs to Type I and the plane is the twinning plane K\. 
If the Miller indices of the 180° rotation axis are rational, the 
twinning mode refers to Type II and the direction of the 
rotation axis is the twinning direction r/^ Since a compound 
twin has two 180° rotations with rational K u K 2 , r\\ and r] 2 , the 
plane normal to the 180° rotation axis that offers the minimum 
shear should be the twinning plane 

In contrast to the SEM/EBSD examination, the TEM 
determination process involves examining the spot diffraction 




Figure 1 

Lattice plane K l with basis vectors U] and u 2 and reduced basis vectors e[ 
and e 2 . 



image (Nishida et al., 2008). For Type I and compound twins, 
the diffraction image - obtained on condition that the incident 
beam is parallel to the Ki plane - consists of two sets of 
reflections that are in mirror symmetry to each other with 
respect to the K { reflection. Thus, the K x plane can be iden- 
tified. For Type II twins, the diffraction image - obtained with 
the incident beam along the i] 1 direction - contains a single 
visible pattern, i.e. the reflections from two twin-related 
crystals overlap each other. The rj 1 direction could also be 
determined. 

Based on the above experimental identification, the other 
twinning elements to define a twinning mode can be further 
derived with the method outlined below. 

2.2. Determination of twinning elements 

2.2.1. Type I and compound twins. According to the clas- 
sical definition, a Type I or compound twin is related to its 
parent by a reflection across the twinning plane K x , where the 
Ki plane is a rational lattice plane with relatively small Miller 
indices. With this condition as starting point, the possible 
twinning direction rji and the magnitude of twinning shear y 
can be deduced in conformity with the minimum shear 
criterion, i.e. the twinning shear that moves all parent lattice 
points to their correct twin positions appears to be the smallest 
in magnitude. Hereafter, our calculations are conducted in the 
direct primitive lattice of the parent crystal. For the coordinate 
transformations between the primitive lattice basis and the 
conventional Bravais lattice basis, we refer to International 
Tables for Crystallography (Hahn, 1996). 

At first, let us choose two basis vectors iii and u 2 in the 
twinning plane K\ and transform them into the reduced 
vectors e t and e 2 , as shown schematically in Fig. 1. The 
reduced basis vectors e t and e 2 must be the two shortest 
translations and the most orthogonal to each other among all 
possible basis vectors in the plane K v Note that such a 
reduced basis is useful for determining the nearest lattice 
point(s) to a given point (not necessarily lattice site) in the 
plane K l . The procedures to find the basis vectors u : and u 2 
and to reduce them to ei and e 2 are detailed in Appendix A 
and Appendix C, respectively. 

Now, we show how to determine the twinning shear vector t 
by use of the reduced basis e 1 and e 2 . Let Plane 0 represent the 
twinning (invariant) plane K 1 that separates the twin lattice 
(above Plane 0) from that of the parent (below Plane 0), as 
shown schematically in Fig. 2. Since the nearest neighbor 
plane (Plane —1) of the parent lattice and its counterpart 
(Plane 1) for the twin lattice are parallel and in mirror 
symmetry with respect to the invariant plane K\, the 
perpendicular projection of Plane —1 onto Plane 1 allows us to 
identify the possible twinning shear vector. Here, we select a 
parent lattice vector OA that ends at the lattice point A on 
Plane —1, and denote by A' the endpoint of the projection of 
vector OA on Plane 1. Obviously, the vector t that joins A' - a 
twin lattice point - to its nearest parent lattice point N on 
Plane 1 defines the twinning direction r\ x and ensures the 
smallest magnitude of shear. The procedures for determining 
the vectors OA and t are described in Appendix B. 
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Furthermore, the interplanar spacing of the twinning plane 
Ki can be easily calculated by the scalar product of OA and m: 

d Ki =\\OX-m\, (1) 

where m denotes the unit vector in the direction normal to the 
twinning plane K x . Thus, the magnitude of shear is given by 

Y=W/d Kl . (2) 

Once the shear vector t and the magnitude of shear y are 
determined, the other twinning elements (rj 2 , K 2 and P) can be 
readily calculated according to the Bilby-Crocker theory 
(Bilby & Crocker, 1965). 

Let I be the unit vector in the twinning direction r\ x and g M a 
vector in the conjugate twinning direction r/ 2 , with reference to 
the parent lattice basis. Applying the twinning operation by a 
shear y along rji, g M is transformed into gM, as shown sche- 
matically in Fig. 3. Since r] 2 is defined by a rotated but 
undistorted lattice line of the shear, g M has the same indices 
(and hence the same length) as g M , if it is referred to the twin 
lattice basis. Moreover, g M and g M lying in the shear plane P 
(perpendicular to K{) are in mirror symmetry with respect to 
the plane that contains the vector V (= d K m) and is perpen- 
dicular to r\\. Thus, the three vectors g M , g M and g form an 
isosceles triangle. As g (= d Ki yl) in the shear direction is 
divided into two equal lengths by V, we obtain 

gM = V-|g = ^ l (m-|yI). (3) 

Notably, g M is not necessarily a lattice vector, and its 
components - expressed in terms of the parent lattice basis - 
can always be transformed into rational indices. Once the 
lattice vector in the r\ 2 direction is determined from g M , the 
shear plane P and the conjugate twinning plane K 2 can be 
easily calculated by the vector cross product (Bilby & Crocker, 
1965). 




OA m 



Figure 2 

Illustration of the nearest neighbor plane (Plane —1) of the parent lattice 
and its counterpart (Plane 1) for the twin lattice, which are parallel and in 
mirror symmetry with respect to the invariant plane K t (Plane 0). The 
twinning shear vector t is represented by the displacement from a parent 
lattice point N to the nearest twin lattice position A' on Plane 1. 



Table 1 

Reciprocal relationship of twinning elements in dual spaces. 
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2.2.2. Type II twin. By definition, a Type II twin is related to 
its parent by a 180° rotation about the twinning direction r] 1 or 
a reflection across the plane normal to the twinning direction 
r}\. Let us first recall the fundamental relationships between 
direct lattice and reciprocal lattice. Every lattice vector in the 
direct space corresponds to a set of lattice planes normal to 
this vector in the reciprocal space, and vice versa. Thus, the 
twin relationship of a Type II twin in the direct space can be 
equivalently expressed by a reflection with respect to the 
plane that is normal to r\i in the reciprocal space, or in other 
words, a Type II twin in the direct space is visualized as a Type 
I twin in the reciprocal space. As the two spaces are strictly 
linked to each other, we can see that, when the direct lattice 
undergoes twinning, the reciprocal lattice is subject to the 
same deformation (shear in the same direction and with the 
same magnitude) and verse visa. In this context, the deter- 
mination of the twinning elements of Type II twins can follow 
the same procedure as that of Type I, except that all the 
calculations should be conducted in the reciprocal space. 
Moreover, the resultant directions (planes) in the reciprocal 
space correspond to the same indexed planes (directions) in 
the direct space, as summarized in Table 1. 

3. Conclusions 

As a widely observed and intrinsic process, crystal twinning 
has a broad impact on the microstructures and properties of 
crystalline materials. So far, the classical theory of twinning 
has advanced greatly the study of twining, but it often suffers 
from insufficient information for practical determination of 
full twinning elements. To progress beyond this state, a general 
method is elaborated based on the minimum shear criterion, 
using the experimentally identified possible twinning plane 



g 




Figure 3 

Transformation of vector g M (in the direction ;j 2 ) into vector gj^ by a 
magnitude of shear y along the direction r)\. Note that g M and gjuhave the 
same length and are in mirror symmetry with respect to the plane 
perpendicular to the K i plane and the shear plane P. 
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for Type I twins or the twinning direction r] 1 for Type II twins 
and the crystal structure as input. As a first step, it determines 
a reduced basis of the invariant lattice plane that serves as the 
mirror plane (in the direct space for Type I twins and in the 
reciprocal space for Type II twins) between the parent and 
twin lattices. Then, a lattice vector - with its origin at the 
invariant lattice plane and its end at the nearest neighbor 
lattice plane of the same set - is selected from the parent 
lattice and projected onto the counterpart lattice plane of the 
twin lattice. Among the vectors that join the endpoint of the 
projected lattice vector to the surrounding parent lattice 
points forming the reduced basis, the shortest vector defines 
the twinning direction and the twinning shear. Finally, the 
other twinning elements can be easily calculated using the 
vector product operations. The present method, as it stands, is 
highly significant for facilitating the study of twinning in a 
variety of crystalline materials. 

APPENDIX A 

Determination of base vectors u, and u 2 on a lattice 
plane 

In crystallography, a lattice plane P with a given Bravais lattice 
is usually described by the Miller indices (hkl), i.e. a set of 
three integers with the greatest common divisor gcd(/i, k, I) = 
±1. Assume |gcd(/z, k)\ = d; then gcd(d, /) = ±1. If an arbi- 
trary lattice vector u with the Miller indices [uvw] lies in the 
plane P, it has 

hu + kv + Iw = 0 or hu + kv — —Iw, (4) 

where u, v and w are integers. Since |gcd(/j,/c)| — d, the 
following relation holds: 

(h/d)u + (k/d)v = -l(w/d). (5) 

Let h/d = h', k/d = k' and w/d = w'\ then ti , k! and w' are 
also integers. Equation (5) can be written as 

h'u + k'v — —Iw' . (6) 

As gcd(h' , k') — ±1, one can find two integers u 0 and v 0 that 
satisfy the following relation according to Bezout's theorem: 

h'u 0 + k'v 0 — ±1. (7) 

Multiplying both sides of equation (7) by (=p/w'), we obtain 

h'(^lw')u 0 + k'(^flw')v 0 = -Iw'. (8) 

Let (^flw')u 0 — u and (=pZw')v 0 = v; then equation (8) becomes 

h'u + k'v = -Iw'. (9) 

Subtracting equation (9) from equation (6), we have 

h'(u -u) + k'(v - v) = 0. (10) 

Since gcd(h' , k') = ±1, there exists an integer a such that 

U — u" = —ak', v — v" = ah' . (11) 

Rearranging equation (11), we obtain 



u = u — kk' = =p/w'w 0 — ak/d, 

v = v + Xh' = =p/w'v 0 + ah/d, (12) 

w = w'd. 

By definition, the basis vectors are a set of linearly inde- 
pendent vectors such that each vector in the space is a linear 
combination of the vectors from the set. Therefore, equation 
(12) proves that the vector (=plw'u Q — ak/d, =F/w'v 0 + ah/d, 
w'd) constitutes the basis vectors of the plane (hkl). Setting 
a = 1 and W = 0, and a = 0 and w' = 1, respectively, we 
obtain two basis vectors: 

Ul = (-k/d, h/d, 0), u 2 = (=fZk 0 , T'v 0 , d), (13) 

where u 0 and v 0 are the Bezout coefficients of equation (7). 
With the Euclidean algorithm, u 0 and v 0 can be easily calculated. 

APPENDIX B 

Determination of lattice vector OA and shear vector t 

B1 . Lattice vector OA 

According to the fundamental law of the reciprocal lattice 
(Authier, 2001), for an arbitrary vector OA with its origin O at 
the zeroth plane of a family of lattice planes (hkl), if it inter- 
sects the wth plane at the point with coordinates (x, y, z), the 
following relation holds: 

hx + ky + lz — n. (14) 

Let OA be the lattice vector with the Miller indices [—2m 
— 2v — 2w] and K l the invariant plane with the Miller indices 
(hkl), as shown in Fig. 2. Then, we have 

— 2uh — 2vk — 2wl = — 2 or uh + vk + wl = 1 . (15) 

The Bezout coefficients u 0 , v 0 and w 0 of equation (15) can be 
calculated with the Euclidean algorithm, and hence the lattice 
vector OA. 

B2. Shear vector t 

Consider a lattice vector OA with its origin at O on Plane 1 
and its end at A on Plane —1, as shown in Fig. 4. 




Figure 4 

Perpendicular projection A' of a lattice point A of Plane —1 onto Plane 1. 
The shortest vector joining A' to the nearest lattice point is denoted as the 
shear vector. 
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Let A' be the perpendicular projection of the lattice point A 
on Plane 1. The shear vector t is defined as the shortest vector 
among all vectors that connect A' with the surrounding lattice 
points on Plane 1. Introducing the reduced basis ei and e 2 , we 
can derive from Fig. 4 that 

OA' = OA+ (OA m), 

O'A = OO' - (Aej + ySe 2 ), 
qA = O'A - e 2 , (16) 
P A = O'A - ( ei + e 2 ), 
rA = O'A - e u 

where m is the unit vector of the plane normal. By comparing 
the lengths of O'A', qA', pA' and rA', the shortest vector t can 
be easily found. 



APPENDIX C 

Transformation of basis vectors u, and u 2 into the 
reduced basis e, and e 2 

To find the closest lattice point to the projection A' and thus 
the minimum shear of the twinning, it is essential to establish a 
reduced basis, i.e. the two shortest lattice vectors that are most 
orthogonal to each other (Zuo et ah, 1995). With the basis 
vectors Ui and u 2 determined according to Appendix A as 
input, the reduced basis e : and e 2 can be derived using an 
iterative procedure, as described below. 

Let e l be the shorter vector between the two base vectors, 
i.e. |e x | < |e 2 |. Then, the new base vectors are derived from 

e'j = e 1? e 2 = e 2 — se l . (17) 

To render the two vectors orthogonal to each other, this yields 

£ = ^X (18) 
e i - e i 



If e < 0.5, e[ and e 2 deliver the reduced basis vectors. 
Otherwise, e is rounded into the nearest integer and the above 
procedure is repeated until e < 0.5. 
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